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1. Introduction 

So far, the only well established hadron states are mesons (qq) and baryons (qqq). In spite of all 
the efforts, no exotic states like tetraquark (qqqq), pentaquark (qqqqq), or hybrid (qqG) states have 
been confirmed beyond doubt. Perhaps the most prominent tetraquark candidate is the Z + (4430) 
resonance, discovered by Belle [1]: it decays to 7T + y/, so it must have a minimal quark content 
ducc, but it has not been confirmed by Babar (yet) [2]. 

Here we address the question whether the lightest scalar mesons a (I = 0), K (I = 1/2) and 
£20(980) (/ = 1) are tetraquarks 1 , as many scientists strongly believe since Jaffe's proposal in 1977 
[3] (see for example [3]). The a resonance (ni a ~ 400 — 600 MeV) is now widely accepted since its 
pole was determined in a model-independent way [5]. The pole for K resonance (m K ~ 600 — 800 
MeV) was determined in a similar manner [6]. Both resonances have been recently experimentally 
confirmed [7], but they remain slightly controversial. 

The observed ordering m K < m ao ^ 9SQ ^ can not be reconciled with the conventional us and ud 
states since m as > m Uc i is expected due to m s > m c \. In the tetraquark case, the 1 = 1 state [ds] [us] 
with additional valence pair ss is naturally heavier than the 1=1/2 state [sJjffifw] and the resem- 
blance with the observed ordering favors the tetraquark interpretation. 

Light scalar tetraquarks have been extensively studied in phenomenological models [4, 3], but 
there have been only few lattice simulations [8, 9, 10]. The main obstacle for identifying possible 
tetraquarks on the lattice is the presence of the scattering contributions in the correlators. The 
strongest claim for a as tetraquark was obtained for m n ~ 180 — 300 MeV by analyzing a single 
correlator using the sequential empirical Bayes method [8]. This result needs confirmation using a 
different method (for example the variational method used here) before one can claim the existence 
of light tetraquarks on the lattice with some confidence. 

We determine a spectrum of states with J PC = ++ , p = and 1 = 0, 2, 1/2, 3/2 on the 
lattice using the variational method with a number of tetraquark sources and sinks. The two- 
pseudoscalar scattering states 7t(k)n(—k), which have back-to-back three-momentum k = N^f-, 
unavoidably couple to our sources with / = 0, 2 and appear in the resulting spectrum. Similarly, 
K(k)n(—k) appear in the 1=1/2, 3/2 spectrums. We don't consider the more challenging 1 = 1 
channel, since there are two towers K(k)K(—k) and n{k)r\(—k) of scattering states. In the non- 
interacting limit, the scattering states have energies \J m 2 + k 2 + J m 2 K K + k 2 . Our main question 
is whether there is any light state with / = in addition to the n(k)n(—k) tower. We indeed find an 
indication for an additional light state in / = channel, which could possibly be interpreted as a 
with strong tetraquark component. We also find an additional light state in / = 1/2 channel on top 
of the K(k)n(—k) tower, and this state could be interpreted as K with strong tetraquark component. 
We also study the exotic repulsive channels with 1 = 2 and 3/2, where no resonance is observed. 
Indeed, in this case we only find the 7i(k)n(—k) and K(k)7i(—k) towers, but no additional light 
states. 



By tetraquark, we have in mind a state that has a dominant tetraquark Fock component. 
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2. Lattice simulation 

We present results from two simulations: 

• An Nf = 2 dynamical simulation with Chirally Improved quarks, three pion masses m n =318, 
469, 526 MeV ( 2 ), a ~ 0. 15 fm and V = 16 3 x 32 [1 1]. "Narrow" Jacobi smearing is applied 
to all quark sources and sinks [11]. 

• A quenched simulation with valence overlap fermions, m % ~ 230 — 630 MeV, a ~ 0.200(3) 
fm and V = 16 3 x 28 [12]. In this case, quarks have point-sources and point-sinks and we 
analyze 300 configurations. 

We use r = 5 tetraquark interpolators in / = and 1 = 1/2 channels: 

Gx = PP' , 2 = £ ViVj , ^3 = £ AiA' t , 4 = [qiCy 5 q 2 ] [qiCy 5 q A ] , 04 = [qiCq 2 ] [qiCq A ] 
i=l,2,3 i=1.2,3 

(2.1) 

where [q^Cy^q^] and [/73G74] are (pseudo)scalar di-quarks and P = qijsq2, Vj = <7i7;<?2, M = 
qi7ij5q2 are conventional cun - ents. Appropriate quark flavor combinations 3 are taken to get 1 = 
and / = 1 /2. 

We use r = 3 interpolators for / = 2 and 1 = 3/2 

0,=PP', 2 = £ ViVl , 3 = £ (2.2) 

1=1,2,3 1=1,2,3 

with flavor content dudu for / = 2 and sudu for / = 3/2. 
We compute the rxr correlation matrix 

Qj(t) = £^(0|^(^0^(0,0)|0)^o = £z, (,l) zj' l) * e - £< " ) ' (2.3) 

x n 

where Z. = (0|^ ; |?i). Like in all previous tetraquark simulations, we omit the single and double 
disconnected quark contractions for / = 0, 1/2. This approximation discards qqqq <-> qq <-> vac 
mixing. Since we are interested in an intermediate tetraquark state with four valence quarks, there 
is even a good excuse to use this approximation in these pioneering studies. 

The extraction of the energies from the correlation functions (2.3) using a multi-exponential 
fit is unstable. Instead, we solve the generalized eigenvalue problem [13] 

C(t)uW(t) = XW(t,to)C(to)d n) (t), (2.4) 
which gives us the energies via X ^ n \i) and the ratios of z\ via u^ n \t) as follows: 



\I#c jV (t)q(t)\ 



,-E^(t-t ) \LkQk(t)u n k (t)\ z) 



(2.5) 



These relations are formally valid at large t, infinite lattice temporal extent T and to < t < 2to [13]. 



2 We analyzed 200 configurations for lighter two pion masses and 100 configurations for the heaviest pion mass. 

3 The current-current interpolators have flavor structure Iduud — uudd + 5 uuuu + jdddd for / = and Eo=« d s sqqu 
for 1=1/2. The prime in Eqs. (2.1,2.2) indicates that two currents may have different flavor structure. The diquark 
anti-diquark interpolators have flavor structure [Ud] [ud] for / = and [sd] [du] for / = 1/2. 
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3. Results for / = and / = 2 



The typical effective masses and the ratios |Z ; - /Zj- | for 1 = 0, 2 are given in Fig. 1. The 
lines display the three lowest energies of %{k)n{— k) in the non-interacting case. 



1=0 m =318 MeV , 5x5 matrix, dynamical simulation 



J—2 m =469 MeV , 3x3 matrix, dynamical simulation 
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Figure 1: Typical effective masses of the eigenvalues A ("'(f) for / = 0, 2 and « = 1, 2, 3. Corresponding 
ratios |Z ( . /Zy | at given « are also shown (J is the largest component). Results for specific quark masses and 
to = 1 in the dynamical simulation are shown. The lines present the energies of non-interacting n(k)n(—k) 
with k = N^ andiV = 0,l,\/2. 
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Figure 2: The resulting spectrum E^ n ' for / = and / = 2 in the dynamical (left) and the quenched (right) 
simulations. Note that there are two states (black and red) close to each other in I = case. The lines present 
the energies of non-interacting n(k)n(—k) with k = N^f- and N = 0, 1, V2. 

In / = case, we find one state with energy close to 7r(0)7r(0), another state with energy 
close to k(^)k{—^) and we also find an additional light state (close to the lowest state). This 
applies for all quark masses and for dynamical as well as the quenched simulation, as shown in the 
resulting spectrum on Fig. 2. Both figures display results for ?o = 1 and for the case when the full 
5x5 matrix (2.1) was diagonalized. We also used to £ [2,4] and diagonalization of all possible 
4x4 and 3x3 sub-matrices and we find that extracted and \z\ /zj | are almost independent 
on these choices. The ratios \Zf >X ' /Zy >1 ' | for excited states depend only mildly on these choices. 
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In / = 2 case, we find one state with energy close to 7r(0)7r(0), another state with energy close 
to n(^)n{—^) and no additional light state (see Figs. 1 and 2). Again, this applies for all quark 
masses, both simulations and for the range of to G [1 , 4] . In this case we use only a 3 x 3 matrix (2.2), 
which is probably not large enough to capture the energy of the k(^-)tc(—^-) exactly (it naturally 
comes out to high). We point out that our intention was not to capture energy of n(^)n{—^) 
correctly, but to verify that there is no light state in addition to 7r(0)7r(0) in the 1 = 2 channel. 

The energies E^ n > in Figs. 2 and 4 were fitted from the eigenvalues X^ n \t) using 4 

X(")(t) = m» [e- E[ " ]t + e- E(n) ^] +A« [ e -E h t e -E P2 iT-t) +e - Ep2 t e - Epi (T-t^ (3 1} 

which takes into account the effect of finite T and periodic boundary conditions for a scattering 
state P\P% = Ktc (see Appendix A of [9]). The time-dependence of the ground state X^'(t) is 
in good agreement with the expectation (3.1) and the fit typically renders A^ n > ~ w^"\ The effective 
masses, plotted in Figs. 1 and 3, are also defined by taking into account the time-dependence (3.1) 
[9]. 



4. Results for / = 1/2 and / = 3/2 

The effective masses and resulting spectrum for / = 1/2, 3/2 are shown in Figs. 3 and 4. The 
conclusions regarding / = 1/2 is completely analogous to the 1 = case above: there is a light state 
in addition to ^(0)7r(0) and K(^-)n(—^-). Results for the exotic 1 = 3/2 channel are analogous 
to results for 1 = 2: there is no light state in addition to K(0)7c(0) and K(^)n(-^-). 



1=1/2 m =318 MeV , 5x5 matrix, dynamical simulation 



1=3/2 111^=469 MeV , 3x3 matrix, dynamical simulation 
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Figure 3: Typical effective masses of the eigenvalues X^ n \t) for / 1/2, 3/2 and n = 1, 2, 3. Corre- 
sponding ratios |Z ( - /Zy | at given « are also shown (J is the largest component). Results for specific quark 
masses and fo = 1 in the dynamical simulation are shown. The lines present the energies of non-interacting 
K(k)n{-k) with k = N^ and N = 0, 1 , s/2. 

We verified again that these conclusions are independent on the choice of 4 x 4 or 3 x 3 sub- 
matrices of the full 5x5 matrix (for 1=1/2 case) and they are independent on the choice of 
foe [1,4]. 

4 The three-parameter fit of the second (first) excited state for / = (/ = 2) was unstable, so we fixed A" = in this 
case. 
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Figure 4: The resulting spectrum for / = 1/2, 3/2 in the dynamical and the quenched simulation. Note 
that there are two states (black and red) close to each other in / = 1/2 case. The lines present the energies of 
non-interacting K(k)lt(-k) with k = N^f- and N = 0, 1,V5- 



5. Conclusions and outlook 

We determined the energy spectrum of states with J PC = ++ and 1 = 0, 1/2, 3/2, 2 using the 
generalized eigenvalue method with a number of tetraquark interpolators at the source and the sink. 
Our simulation is done at several values of m n and at a single lattice size L. Our main question is 
whether we find any light states in addition to the towers of n(k)n(k) or K(k)n(k) (k = 0, ..) 
scattering states. 

We do find additional light states in the 1 = and 1=1/2 channels. We conclude that these 
additional states may be possibly related to the observed resonances a and K. Since we used 
only the tetraquark sources and discarded the annihilation diagrams, our simulation also gives an 
indication that a and K may have a strong tetraquark component. 

We point out that results from our simulation would have to be confirmed by an independent 
lattice simulation, before making firm conclusions. At this point, we can not completely exclude the 
possibility that the additional light states observed in / = 0, 1/2 channels are some kind of unknown 
artifact from the generalized eigenvalue method, omission of annihilation diagrams, etc. However, 
our results behave well in a number of extensive checks we made to exclude this possibility. 

The ultimate method to study a and K on the lattice would involve the study of the energy 
spectrum as a function of the lattice size L. The resonance mass and width could then be determined 
using the Liischer's finite volume method [14]. The importance of the tetraquark component could 
be determined by comparing {0\^^\n) and (0\qq\n). 
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